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Abstract 

An attempt is made in QCD to explain the growth of total cross-sections 
with energy, without violating the Froissart bound. This is achieved by com- 
puting the phase shifts of elastic scatterings of partons rather than their am- 
plitudes. To render that possible a general formalism of phase-shifts in QCD 
is developed. Computed to two-loop order, agreements with hadronic and 
virtual-photon total cross-sections are quite satisfactory. Predictions for the 
slower rate of growth at higher energies are also presented. 
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The recent observation of the rapid rise of parton density at small x has generated 
much theoretical interest ||. This rise is equivalent to an increase of the total photon-proton 
cross-section 0t(7*p) with its cm. energy y/s. For photons with virtuality Q 2 in the range 
of the HERA data, (Jt{i*p) ~ s a with the power a substantially larger than the power 0.08 
Q observed in hadronic and real-photon-hadron total cross-sections, but substantially less 
than the power ~ 0.5 given by the BFKL Pomeron ||. See Fig. 2. A power growth violates 
the Froissart bound ot ~ (Ins) 2 so eventually all these increases have to slow down. In 
terms of parton density, this means saturation at small x when the partons begin to overlap 
in transverse dimensions 0. 

A considerable amount of work exists in computing parton densities at small x for large 
nuclei, and in understanding its saturation f?|,|8|] . We shall approach this problem by look- 
ing at the total cross-section 0r(7*p), which at least for Q 3> ^-qcd is also amenable to 
perturbative QCD calculations. At sufficiently high energies, the energy variation of total 
cross sections is expected to be universal, so for simplicity we shall only study the energy 
dependence of quark-quark interactions. 

Optical theorem relates total cross-section to the imaginary part of forward elastic scat- 
tering amplitude. If the latter is computed via the exchange of a BFKL Pomeron [||, then 
<7t(s) violates the Froissart bound f|]||. Multi-reggeon exchanges are necessary to restore 
unitarity but that calculation is very difficult |J. We take a different approach by looking 
at the phase shift, which will almost certainly guarantee the Froissart bound. The idea of 



using phase shift is not new |10|, though it remains a challenge to calculate it within the 
framework of QCD. We shall discuss how this can be carried out, then use it to compute the 
quark-quark scattering phase-shift to two-loop order, in the leading-log approximation. The 
result compares well with the energy increase of the HERA data at different Q 2 (Fig. 2), 
and it also agrees with the hadronic data up to laboratory energy of 10 9 GeV (Fig. 3). 

Let A(s, A) be the elastic scattering amplitude at momentum transfer A. At high ener- 
gies it is well described by the impact-parameter representation 
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AO, A) = is J d 2 be iKrb (l - e 2 ^ s ' 6 >) , (1) 

where A and b are two-dimensional vectors in the transverse plane. At large impact parame- 
ters the phase shift 5(s, b) is small and it goes approximately like a exp(-fib) if the interaction 
range is Only for b less than the effective radius R(s) does the phase shift contribute 
substantially to the amplitude, so R can be estimated from the condition S(s,R) ~ 1. If 
a is an increasing function of s, then R(s) and hence ot ~ ttR(s) 2 also increase with s. 
In particular, if a ~ s for I > 0, then R(s) ~ (£/fi)lns, and oy ~ (Ins) 2 reaches the 
Froissart bound. This conclusion is stable and robust, and is qualitatively independent of 
the magnitude of I and \i. Unless a grows faster than a power of s, the Froissart bound is 
almost certainly guaranteed when computed via the phase shift. 

Feynman diagrams tell us how to calculate perturbative amplitudes but not directly 
the phase shifts. This is because the impact-space amplitude A(s, b) = e 2tS ^ s ^ — 1 = 
— X)^Li [2id>(s, b)] n /n\ is given by an infinite sum of all powers of the phase shift. Even if S(s, b) 
is computed to the lowest order, the resulting A(s, A) already contains terms of all orders. 
Thus a proper understanding of phase shifts cannot be obtained unless an infinite sum of 
Feynman diagrams is considered. More generally, if the phase shift is expanded in powers of 
the coupling constant g 2 , 5(s,b) ~ J2m g 2m 5^ m \ then from ([!]) the (2n)th order contribution 
to A(s, b) is given by a sum of products of the phase shifts S^ m \ Individual Feynman diagrams 
certainly do not factorize in this manner so even at a given order we have to sum over many 
Feynman diagrams. A hint for what to sum is given by the formula 5(s, b) ~ / V(s, b, z)dz 
for potential scattering, where 5(s, b) is given by an accumulation of interactions (vertices) 
along the path without involving any energy denominators or propagators. Recall that in 
perturbation theory the energy denominator 1/AE comes from the uncertainty relation, 
and the fact that the time At ~ 1/AE allowed for each interaction in a given Feynman 
diagram is constrained by the times of its two neighbouring interactions. To get rid of 
the energy denominator we must relax the time constraint by summing over diagrams with 
all possible orderings of the interactions. When this is carried out with the help of the 
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eikonal formula [TT| for potential scattering and for QED, indeed the energy denominators 
disappear, factorization occurs, and the phase shift is given by the first Born approximation. 

For QCD, the presence of non-commuting colour matrices makes factorization consider- 
ably more difficult to achieve, and representations in terms of phase shifts harder to obtain. 
In what follows, we shall first discuss how this is done for the simpler case of scattering 
from an external colour source, before tackling the more difficult problem of quark-quark 
scattering where the leading-log approximation has to be invoked. 

The desired factorization for a fast parton of momentum p interacting with external 



colour sources, in the tree approximation, comes from the decomposition theorem [12.13 



k— 1 k 

A n (b) = £ c mi (b)c m2 (b) ■ ■■e mk (b)/ nE% (2) 

{m} 1=1 j=i 

with the first sum taken over all rrii > 1 so that J2i=i mi — n, and over all k. For n > 1, 
the irreducible amplitude C n (b) is identical to the Feynman amplitude A n {b), except the 
product of colour matrices A ai A a2 ■ • • A an _ 1 A an in A n is replaced by its nested commutators 
[A ai , [A 02 , [• ■ ■ , [A an _ i; A a J]]] in C n . For n = 1, by definition C\ = A- Note that A n is given 
by a product of colour matrices so it contains many colours, but C n is given by their nested 
commutators so it carries colour only in the adjoint representation, same as the colour of a 
single gluon. 

With this factorization it is possible to sum up A n (b) to obtain an impact-parameter 



representation [13], from which the formula for the phase shift can be extracted. The phase 



shift defined in ([]]), with s replaced by 2p, is given by 

where the ellipses consists of commutators of C m of order g 10 and above. The 0(g 10 ) expres- 



sion is explicitly known [|T^] and the higher order expressions can be computed when they 
are needed. 

For quark-quark scattering the presence of loops and colour matrices of the other quark 
makes the problem so much more difficult that we find it necessary to invoke the leading- 
log approximation. It is cruicial to recognize however that the leading-log approximation 



used here may not be the same as what is found in the literature. For example, it is 
known that the BFKL Pomeron |J computed in the leading-log approximation violates 
the Froissart bound; subleading-logs supplied by multiple-reggeon exchanges are needed to 
restore unitarity. This however refers to leading and subleading logs with respect to a fixed t- 
channel colour (singlet in this case). The leading- log approximation to be used in the phase 
shift calculation keeps leading-log terms with respect to fixed colour structures. Leading 
logs in fixed colour structures may give rise to subleading logs in fixed t-channel colours, 
like those supplied by multiple reggeons, which is why leading-log for fixed colour structures 
may give rise to phase shifts that preserves the Froissart bound. Indeed, there is a very close 
relationship between colour structures and reggeon structures that will be discussed further 
below. 

The colour structure of a diagram is defined by its colour factor when quarks carry 
arbitrary SU(N C ) colours. Colour factors are linearly dependent only if they can be related 
through SU(N C ) commutation relations alone, with coefficients independent of N c and the 
specific quark colours. The six diagrams in Fig. 1 have linearly independent colour structures. 

The colour structure of a diagram is said to be primitive if all the lines in the diagram 
remain connected after the top and the bottom quark lines are removed. Using G to denote 
the colour matrices of a colour structure, G a and Gd in Fig. 1 are primitive and all the others 
are not. 

Decomposing gluons attached to one of the quark lines using (||), it can be shown in the 
leading-log approximation [0,0 that (i). Non-primitive colour factors may be considered 



as commutative products of primitive colour factors. For example, using G to denote the 
colour structure of a diagram, then Gb = G 2 a , G c = G a , G e = Gf = G d ; (ii). This colour 
factorization also leads to an amplitude factorization in the impact-parameter space. If the 
impact-parameter amplitude for a primitive colour structure Gi is dj, then the spacetime 
amplitude for the colour factor M(Gi), where M is a monomial function of the GVs, can be 
factorized into M(efj); (iii). As a result of this factorization, diagrams of all orders can be 
summed up to an impact-parameter representation with a phase shift given by 25 (s, b) = 



J2idiGi, where the sum is taken over all linearly independent primitive colour factors Gj. 
To two-loop order in the leading-log approximation, the quark-quark scattering ampli- 



tude A(s, A) is known to be fl6| , [T7|j 



A{s, A)/s = {-g 2 h(A) + (g 4 N c lns/4n)I 2 (A) 
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\\(2ny(k 2 ± + ^), 
S 2 (a-J2 ki±) = I d 2 b e iA * [K ^b)/2n] n . 



i=l 



(4) 



The infrared cutoff \x is introduced as a parameter to simulate hadronic size and confinement 
]TS[. The factorization property mentioned above can be explicitly verified in eq. (0). The 
resulting phase shift to two- loop order is then given by 

25(s,b) = 25 a (s,b)G a + 25 d (s,b)G d , 
26 a (s,b) = -g 2 K (x)/27r + g 4 N c \nsK 2 (x)/lQ7i 3 

-g 6 N 2 \n 2 sV(x)/32tt 2 , 
25 d (s, b) = ig 6 In s [k^(x)/1Qtt 4 - V(x)/4tt] , (x = fib) 
A 2 1 2 (A) = J d 2 b e iK - $ V(fib). (5) 

The colour factors are respectively G a = \ a x A a and G d = A a Ab x \ c ^d(iface)(ifbed), where 
x indicates tensor product of colour matrices associated with the two quark lines. To 
obtain the elastic amplitude and hence the total cross section, we must extract from ([!]) 
the contribution of the colour-singlet. The resulting formula for the total cross section at 
N r = 3 is 
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a T (s) = ^ 2 Jd\[l- 2 -e- 



cos 





COS 



( 



46, 
3 



a 



) 



(6) 



where the damping exponent D = —i5d is real and so is 5 a . Note that these two functions 
depend on the impact parameter only through the combination x = fib, so that all the \x 
dependences are factored out into the l//x 2 factor in front of the integral. In this way the 
uncalculable parameter [i affects only the overall magnitude of the cross-section, but not 
its energy dependence which is expected to be universal (as is experimentally the case for 
hadronic and real photon total cross-sections). We may therefore use this formula for quark- 
quark total cross-section on the energy dependence of experimentally accessible beams and 
targets. 

To compare the prediction of @ with experiment we have to realize that perturbative 
QCD contains no energy scale per se, so an energy scale A(Q) must be introduced externally 
to replace all Ins factors by ln(s/A 2 (Q)). Recall that (||) is obtained in the leading-log 
approximation, where Ins comes from an integral of the form ln(s/A 2 ) = f A 2i s dw/uj, with 
A being a combination of the other energy scales present (masses, momentum transfers, 
virtualities) . For massless quarks and forward amplitude, on dimensional grounds it is 
therefore reasonable to assume A(Q) = A +cQ, with A to be of the order of Aqcd, which we 
shall take it to be 0.2 GeV, and a parameter c which must be determined phenomenologically. 
We shall also take the coupling constant to be a s = g 2 /47r = 0.26, which is its value 
at Q 2 = 0.65 (GeV/c) 2 , the largest Q 2 involved in these data. The dependence on this 
parameter is very weak so its actual value is really not that important. In particular, we 
will not let it vary with Q 2 though there is no problem to include such a variation. As 
mentioned before, the parameter fi simulates confinement and hadronic-size effects so it is 
not computable within the framework of perturbative QCD so it must be fitted for each Q. 

Total cross sections for deep-inelastic |1|] and hadronic |]19[ data are shown in Figs. 2 
and 3 respectively. The energy variations predicted by (|5[) and (|6|) are given by the solid 
curves, with c = 4. No attempt has been made to obtain a best fit. Shown for comparison 
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in Fig. 2 are also the dotted lines representing a power variation s 08 and a dashed line 
giving s ' 5 . In Fig. 3 the dashed line represents the Donnachie-Landshoff fit [f|,[19| cr<r = 
22s a079 + 56.1s-°- 46 mb. 

At s values above a few tens of (GeV) 2 for the deep-inelastic data, and perhaps a bit 
higher for the hadronic data, the energy variation of this theory, with essentially one free 
parameter c, agrees well with the data. The s-variation of this theory is slower than a 
power, as it must be to obey the Froissart bound, but at this energy range it is really not 
very different from a power as can be seen in the plots. This theory with its leading-log 
approximation is not designed to work at low energies, but it is amusing to note from Fig. 3 
that the dotted line with the low-energy correction term s -0 - 46 neglected will also come to 
very similar values at s ~ 10 2 (GeV) 2 . 

To summarize, we have suggested that the energy dependence of total cross sections 
should be understood through the phase shifts of the elastic amplitudes, because the Froissart 
bound is almost certainly guaranteed in such an approach. The problem of how to calculate 
the phase shifts in QCD at high energies is discussed, exactly for partons scattered from 
external sources, and in the leading-log approximation for parton-parton scatterings. The 
latter is also explicitly calculated to two-loop order and compared with experiments. With 
an energy scale A(Q) = Aq + cQ this theory compares well with the energy variations of 
both the deep inelastic and the hadronic data. 

This research is supported by the Natural Sciences and Engineering Research Council 
of Canada. I wish to thank Frangois Corriveau for help in obtaining the experimental data 
and Greg Mahlon for instructing me how to plot them. 
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FIGURES 



(a) (b) (c) (d) (e) 



FIG. 1. Colour structure diagrams. 
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FIG. 2. 7*-proton total cross-sections as a function of s. Data are taken from Ref. p] where 
references to the original experiments can be found. Q 2 of the photon, in (GeV) 2 , are given on the 
right. Where two numbers are listed, the second one refers to the low-energy data of E665. The 
solid curve is the prediction of the present theory. The dotted curve and the dashed curve depict 
respectively energy variations of s - 08 and s ' 5 . 



10 



200 




Piab (GeV/c) 



FIG. 3. pp and pp total cross sections as a function of the laboratory momentum pi a b- Data are 



taken from Ref. [19], dashed line is the Donnachie-Landshoff fit, and the solid line is the prediction 
of the present theory. All parameters are identical to those used in Fig. 2. 
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